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ABSTRACT 



Context. We observed the Crab pulsar in October 2008 at the Copemico Telescope in Asiago - Cima Ekar with the optical photon 
counter Aqueye (the Asiago Quantum Eye) which has the best temporal resolution and accuracy ever achieved in the optical domain 
(hundreds of picoseconds). 

I I . Aims. Our goal was to perform a detailed analysis of the optical period and phase drift of the main peak of the Crab pulsar and 

HH ■ compare it with the Jodrell Bank ephemerides. 

Methods. We determined the position of the main peak using the steepest zero of the cross-correlation function between the pulsar 
signal and an accurate optical template. 

Results. The pulsar rotational period and period derivative have been measured with great accuracy using observations covering only 
Q ■ a 2 day time interval. The error on the period is 1.7 ps, limited only by the statistical uncertainty. Both the rotational frequency and its 

5-H ' first derivative are in agreement with those from the Jodrell Bank radio ephemerides archive. We also found evidence of the optical 

j/j , peak leading the radio one by ~ 230 fis. The distribution of phase-residuals of the whole dataset is slightly wider than that of a 

d ■ synthetic signal generated as a sequence of pulses distributed in time with the probability proportional to the pulse shape, such as the 

average count rate and background level are those of the Crab pulsar observed with Aqueye. 

Conclusions. The counting statistics and quality of the data allowed us to determine the pulsar period and period derivative with great 
accuracy in 2 days only. The time of arrival of the optical peak of the Crab pulsar leads the radio one in agreement with what recently 
reported in the literature. The distribution of the phase residuals can be approximated with a Gaussian and is consistent with being 
completely caused by photon noise (for the best data sets). 
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1 . Introduction IZampieri et alj 20Tlb. des pite the secular decrease of the lumi 



. nositv (^asut i et aT 19961) and the presence of glitches and tim 



X 



Since its discovery dStaelin & Reifensteinlll968t IComella et alj ing noise. Occasiona lly small variations of the shape of the pulse 

1 19691) . the pulsar in the Crab nebula has been one of the most have been observed dKarpov et alj|2007l) . 
targeted objects in the sky at all wavelengths, from radio to very Several issues concerning the pulsar engine and the actual 

high energy-rays, serving as a test bed for pulsai- theories as well geometry of the emission regions are still debated, ranging from 

H , as for studying astrophysical non-thermal processes. Opti cal pul- the nature and location of the acceleration mechanism, to wave- 

- - ■ sations were discovere d more than 40 yeai's ago (Cocke et alj length dependent variations of the pulse profile, to properties of 

|1969|:|Lynds etaL||l969i) and the Crab pulsar was indeed the first the Giant Radio Pulses (GRPs). In particular, the study of GRPs 

celestial object to be detected as a pulsating source in the optical is cuiTently a very active field, with deep theoretical implica- 

band. tions. So far GRPs have been observed in a handful of pulsars. 

The optical light curve of the Crab pulsar has been monitored However, an analogous optical phenomenon has been observed 

through the yeai's using a vaiiety of telescopes and instruments only at the Crab. GRPs seem to show a weak correlation with op- 

(e.g. Slowikowska et al. 2009). The pulse shape is characterized tical pulses dKuiper et aL||2003|; |Rots et aL||200l), which are on 

by a double peak profile, separated in phase by ~ 140°. The average 3% brighter when coincident with GRPs (Shearer et alj 

shape is similar through the entire electromagnetic spectrum, 2003 ). Recent coe v al tim ing at optical and radio wavelengths by 

although the morphological details differ substantially from ra- Oosterbroek elalj (|2008|) found a 255 ± 21 /^s delay of radio 

dio to gamma-rays. Wavelength-dependent changes in the pul- with respect to optical pulse. 

sar properties have been reported also by Percival et al. (1993) In the last few years we started a monitoring programme 

(peaks width and separation larger in the V band than in the UV) of the Crab pulsar in the optical band aimed at studying the 

and bv iFordham et alJ(l2002h . The pulse shape is verv stable (e.g. long term stability and sub-yus structure of its pulse shape, and 
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Table 1. Log of October 2008 Crab pulsar observations per- 
formed with Aqueye mounted at the 182cm Copernico telescope 
in Asiago. The start time of the observations is the GPS integer 
second, accurate to ±30 nanoseconds. 



Starting time Duration 
(UTC) (s) 



1 


October 


10, 


23: 


:45: 


:14 


898 


2 


October 


11, 


00: 


:05: 


:07 


1197 


3 


October 


11, 


01 


:00: 


:22 


1797 


4 


October 


11, 


01 


:45: 


:44 


1797 


5 


October 


11, 


02: 


:23: 


:07 


1631 


6 


October 


11, 


03: 


:23: 


:46 


1197 


7 


October 


11, 


23: 


:08: 


:03 


292 


8 


October 


11, 


23: 


:25: 


:09 


3597 


9 


October 


12, 


00: 


:54: 


:31 


1794 


10 


October 


12, 


23: 


:03: 


:59 


292 


11 


October 


12, 


23: 


:13: 


:57 


3998 


12 


October 


13, 


00: 


:57: 


:07 


7194 



at performing accurate optical timing of the main peak. By 
comparing the behaviour in other wavebands, especially radio, 
we aim to improve the understanding of the geometry of the 
acceleration site. Observations are performed by means of a 
very fast single photon-counter instrument, Aque ye, mounted 
at the 182 cm Copernico Telescope in Asiago (B arbieri et al.l 
[2009; Zampieri et al. 2011). The design of Aqueye follows 
that of Quan tEYE (the Quantum EYE; iDravins et alJ l2005t 
iBarbieri et alJ l2008i) . an instrument specifically tailored for 
studying rapid optical variability of astrophysical sources with 
the ESO E-ELT. A second version of the instrument, named 
Iqueye, ha s been installed and successfully used at the ESO NTT 
telescope dNaletto et al.ll2009ll2010j) . 

In a preliminary investigation ("Zampie ri et al.l20TTI) we con- 
centrated on the pulse shape stability of the Crab pulsar and 
found that it is stable at the level of ~ 1% on a timescale of 14 
years. This result reinforces evidence for decadal stability of the 
inclination angle between the spin and magnetic axis, and of the 
thickness of the emission region. Here we present a follow-up 
investigation reporting accurate phase analysis of optical timing 
of the main peak and comparing it wi th radio ephemeri des of the 
Jodrell Bank (JB) radio Observatory (iLvne et al.lll993h . 

The plan of the paper is as follows. In section |2] we list ob- 
servations of the Crab pulsar performed with the Aqueye instru- 
ment and discuss the barycenterization procedure. In Section [3] 
we illustrate the implementation of phase analysis. In Section|4] 
our results are presented and compared with radio ephemerides. 
Section|5]summarizes the conclusions. 



2. Observations 

The Crab pulsar was observed with Aqueye mounted at the 182 
cm Copernico Telescope in Asiago. The observations were per- 
formed in 2008 and lasted for three nights, starting from October 
10. The sky was clear and seeing conditions fair (1.5 arcsec aver- 
age). For a timing log of observations, see Table[T] During each 
observing run we recorded the arrival time of ~ 0.15 - 3.6 x 10^ 
photons, time-tagged with a relative time accuracy of ~100 pi- 
coseconds and an absolute precision (referred to UTC) better 
than 500 ps (for d etails about t he timing accuracy of the acqui- 
sition system see iNaletto et al.ll2b0ft) . To our knowledge, this is 



Table 2. Geocentric coordinates of the 182 cm Copernico tele- 
scope in Asiago. The 3cr uncertainty is 0.3 m. 



X 


y 


z 


(m) 


(m) 


(m) 


4360966.0 


892728.1 


4554543.1 



the most accurate measurement of photon arrival times from the 
Crab pulsar ever obtained in the optical band. 

The time-tag of each detected photon in the unbinned time 
series was reduced to the Solar System barycentric time using 
the software Tempo^] ( Hobbs et al. 2006; Edwards et al. 2006b. 
The adopted position of the Crab pulsar is that reported in the 
Jodrell Bank monthly ephemerides (R.A. 05h 34m 31.97232s, 
DEC. -h22° 00' 52.0690" [J2000]), with no correction for proper 
motion. To perform this conversion the software needs also accu- 
rate value of the observatory geocentric coordinates. They were 
obtained with a GPS receiver, which was connected to an an- 
tenna with a length compensated cable and situated at the dome 
of the telescope. At least 6 GPS satellites were used in the po- 
sitional data acquisition, which typically lasted 3 hours and was 
repeated for several days. Finally the position of the antenna was 
referred to the intersection of the telescope hour angle and dec- 
lination axes by laser assisted metrology. We estimate the posi- 
tional error to be ~30 cm, amply sufficient for the purpose of this 
paper. 

In order to compare our ephemerides to those reported in the 
JB Observatory radio archive we baricentered the time-tags also 
in the Tempo 1 emulation mode. In doing so, we found an error in 
the actual value of the Roemer delay, caused by some inconsin- 
stency in the Earth configuration files related to the calculation 
of the polar motion of the EartlQ. This problem was solved by 
using updated Earth configuration files loaded from the Tempo2 
SCM repositor>Q. 



3. Phase-analysis 

We analyzed the evolution of the phase of the main peak of the 
Crab pulsar using the bary centered event list. A reference period 
Pinit is assumed and the light curve is divided into n seconds long 
segments. Each segment is then folded over P,„„ and is binned 
at ~ 1/300 in phase. The phase of the main peak is determined by 
cross-correlating the pulse shape with a template, as summarized 
in the App endix. Our method is conc eptually similar to the one 
adopted by 'Oosterbroe k et al.l ( l2008l) and is more accurate than 
previous approaches based on fitting the main peak with a simple 
analytic funct ion like a Lor entzian, Gaussian or parabola (see 
e.g. , Oosterbroek et alj|2006l) . given its asymmetric shape. 



' http://www.atnf.csiro.au/research/pulsar/ppta/tempo2 
- The Roemer delay computed in Tempo2 showed an anomalous os- 
cillation at around the time of our observations, reachin g a maximum 
value of 30 /js (instead of i 35 ns; see e.g. [Hobbs et al. 2006). This 
caused a drift of the phase of the main peak of ~1 ms per day and a 
lengthening of the rotational period of ~0.4 ns. 
^ http://tempo2.cvs.sourceforge.net/viewvc/tempo2/tempo2/T2runtime/ 
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Following the standard pulsar spin down model, we describe 
the phase drift of the main peak with respect to uniform rotation 
using a third-order polynomial, i.e.: 

1 9 1 ^ 

A0(f) = 0(0-0 (0 = 0O+(v-V,„;,)(f-fo) + -v(f-fo) +7Kf-fo) (1) 

/ o 

where fo is a reference time, 0o - 4>ik)) is the phase of the main 
peak at fo, (p'{t) - {t - tQ)/Pi„it is the phase for constant rotation 
at frequency v,„„ - 1 /P,„,y, and v, v, v are the actual rotational 
frequency and its first and second derivatives, respectively. P,„,-, 
is chosen in such a way that A(f> varies slowly during the period 
of observation. For a 2 days baseline, the linear and quadratic 
terms, i.e. the first and second derivatives of the phase, are suf- 
ficient to describe the drift, and we can safely neglect the cubic 
term in equation ([1]). Then, the expression for A0(f) (eq. HI) be- 
comes of the form 

0(0 = 00 + a{t - to) + bit - tof , (2) 

where 0o, a (in units of s^^) and b (in units of s^^) are determined 
by fitting 0(0 to the observed phase-drift. After determining 0(0 
from the fit, the phase of the main peak is given by: 

0(0 = 0'(O + A0(O = Vi„i,(t - fo) + if>{t) . (3) 



30000 



25000 



20000 



15000 




0.5 1 1.5 2 

phase 



Fig. 1. Folded light curve of the Crab pulsar as a function of 
phase for the Aqueye observation 4 in Table [T] The folding pe- 
riod and the bin time are 0.0336216417 s and 33.6 ;us, respec- 
tively. The typical double peak profile of the pulse is recogniz- 
able. For sake of clarity two rotations of the neutron star are 
shown. Phase zero/one corresponds to the position of the main 
peak in the radio band and is marked with a vertical dashed line. 



3.1. The radio phase from the Jodrell Bank ephemerides 
archive 

We compared the phase of the Crab pulsar measured by Aqueye 
with that reporte d in the radio archive at the JB ObservatorjQ 
dLvne et al.ll993h . The phases of the JB ephemerides are those of 
the main peak at infinite frequency at the barycenter of the Solar 
System. The observed barycentric radio phase 0r(f) is obtained 
using the values of v^, and Vr nearest to our observing epochs 
and reported in the archiv^. The radio phase drift is given by an 
expression similar to equation ([T]): 

1 7 1 ^ 

A0, = 0r(O-0 (0 = 0,-,o+(vr-v,>„Y)(f-fo)-H-v,(f-fo)^-i--v,(f-fo)' 

/ 6 

The optical phase is in agreement with the radio one if the phase 
drifts inferred from equations ([T]i and ^ are in agreement. The 
radio phase 0r,o at epoch fo is calculated by means of a quadratic 
extrapolation starting from the closest value in time reported in 
the JB radio ephemerides (Oct 15, 2008) and is in agreement, 
within the errors, with the value calculated with a Fortran code 
available at the JB radio ephemerides website. No known glitch 
(down to the intensity to which radio monitoring is sensitive to) 
affected the Crab in the interval of time between the determina- 
tion of the Crab radio parameters and the optical observations. 



~ 5500 counts s^'. The count rate of the background, estimated 
from the ofF-pulse region of the folded light curve (see FigurelTJ, 
is ~ 4500 counts s"', which implies a total number of net source 
photons of ~ 2.6 X 10^, time tagged to better tha n 500 ps with re - 
spect to UTC. As we already point out in Zamp ieri et al.l (1201 lb . 
the pulse shape agr ees w ell with the 33 years-old pulse profile 
ob tained bv | Groth| ( 1975h . as well as with the more recent one 
bv lKarpov et all (I2007h . 



4. 1 . Opticai phase drift 

{4-Jn the following we focus on the detailed analysis of the opti- 
cal phase drift of the main peak of the Crab pulsar and com- 
pare it with the behaviour observed in the radio using the JB 
ephemerides. Figure |2] shows the best-fitting parabola to the 
phase-drift of the Crab pulsar measured by Aqueye. The phase of 
the main peak is calculated using intervals 2 seconds long. The 
typical 1-cr uncertainty on the position of the peak is ~ 30yus. 
The best-fit gives a reduced ~ 1.06. The reference epoch fo 
is MJD=54749.0, while the reference rotational period P,„„ = 
l/vi„i, used to fold the light curve is P,„/, = 0.0336216386529 s. 
The best-fitting parabola of the optical phase drift 0(f) (eq. 121) 



4. Results 

In Figure [1] we show the light curve of the Crab pulsar folded 
over the average spin period for one of the Aqueye observations 
(see also Zampieri et al. 201 1). The light curve includes the con- 
tribution of both the pulsar and the nebular background entering 
the Aqueye pinhole entrance aperture. For a whole observation 
the counting statistics is large and hence the bin time is smaller 
than that adopted for the phase analysis, for which the typical 
integration time is a few seconds (see below). The average count 
rate of the Crab pulsar (all channels) measured by Aqueye is 



0(0 = (1.021431 ±0.000081) 

- [(3.21329 + 0.0001 1) X 10"^r'](f-fo) 

+ [(1.859380 ± 0.000029) x 10"'° s'^] (f - fo)^ . 



(5) 



http://www.jb.man.ac.uk/~pulsar/crab.html 
' http://www.jb.man.ac.uk/~pulsar/crab/all.gro 



The quoted uncertainties are the 1-cr errors for one interesting 
parameter. Figure |2] shows also the phase residuals remaining 
after subtracting equation (|5]l from the measured phase drift. The 
distribution shows a spread of ^ 100 /is (~ 0.003 cycles) and is 
rather symmetric around zero (see Figure HJi, testifying that all 
the observations are phase-connected. 

The rotational frequecy is v = d<p{t)/dt. Thus, inserting equa- 
tion (|5]l into equation (O and taking the derivative with respect 
to f, we obtain an expression for the rotational frequency v and 
period P = 1/v at a given (bary centered time) f. Table [3]lists the 
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1.0 1.5 2.0 2.5 3.0 



time since MJD=54749.0 (days) 

200| , , , , 



150 




-150 



1.0 1.5 2.0 2.5 3.0 

time since MJD=54749.0 (days) 

Fig. 2. Top panel: Phase-drift of the main peak of the Crab pul- 
sar (changed sign) measured during the observing run in Asiago 
in October 2008. The (red) curve is the best-fitting parabola 
(eq|l2l). The reference epoch fo is MJD=54749.0, while the ref- 
erence rotational period is P,„,, = 0.0336216386529 s. Bottom 
panel: Phase residuals (in jUs) after subtracting the best-fitting 
parabola to the phase-drift. 



rotational periods of the Crab pulsar measured with Aqueye and 
those from the JB ephemerides at 3 reference dates (barycen- 
tered MJD=54750, 54751, 54752). The differences between op- 
tical and radio ranges from ~ 2 to 4 ps. The statistical error on 
the optical rotational periods is 1.7 ps (Icr error). The quoted er- 
ror from the JB radio archive is ~ 0.1 ps. Thus, the rotational 
periods calculated by Aqueye agree within the statistical error 
with those derived from radio measurements. Also the measure- 
ments of the first derivative of the rotational frequency v are in 
agreement within the errors: VA^Hcyc - 3.71876x 10"'"±6x 10"'^ 
s-2 and vjB = 3.718655 x lO^'" + 2 x 10""' s^^. 

We emphasize that, in order to compare our data with those 
of the JB radio ephemerides, we barycentered the time-tags of 
the optical photons in Tempo2 using the Tempo/Tempo 1 emu- 
lation mode. This is because, for historical reasons, the JB ra- 
dio ephemerides are calculated using Tempo, which is the older 
version of the software used for barycentering. The systems of 
time adopted in the two packages are different. Tempo2 uses the 



Table 3. Rotational periods of the Crab pulsar measured by Aqueye in 
2008 compared to those reported in the Jodrell Bank radio ephemerides. 
The time- tags were barycentered in Tempo 1 emulation mode. 



MJD" /'(Aqueye)* P (JB) 

(s) (s) 

54750.0 0.033621638649 0.033621638653 

54751.0 0.033621674970 0.033621674973 

54752.0 0.033621711290 0.033621711292 



'' MJD at the solar system barycenter (Tempo 1 mode). 
crp= 1. 7 ps (68% statistical error) 



Barycentric Coordinate Time (TCB), while Tempo the barycen- 
tric dynamical time (TDBfl In Table |4] we report rotational pe- 
riods after barycentering with Tempo2 in TCB units. Although 
there are other differences between the Tempo/Tempo 1 mode 
setup used for calculating the JB ephemerides of the Crab and 
the full Tempo2 mode (as for example the adopted Solar System 
ephemerides), the main difference between rotational periods re- 
ported in Tables[3]and|4]is due to the use of TCB units (SI units) 
instead of TDB. The rotational periods of the Crab are ~ 0.5 ns 
longer than those measured using Tempol in TDB units. If we 
take the ratio of the periods in Table |4] (TCB units) to those in 
the second column of Table |3] (TDB units) we find that the ratio 
of the two time units in K ~ I + 1.53 x 10"^ + 1.3 x 1 0~'°, con- 
sistent with the value reported by Irwin & Fukushima ( 1999) and 
jiobbs et al. (2006). The constant K sums up a contribution from 
the linear term of the Einstein delay, Lc, and another term from 
the gravitational plus spin potential of the Earth, Lg. Thanks to 
its timing capability and performances, Aqueye can put in evi- 
dence the occurrence of the corrections Lc and Lc to the pulsar 
spin period in only two days of data taking on the Crab pulsar. 



Table 4. Rotational periods of the Crab pulsar measured by Aqueye 
in October 2008. The time- tags were barycentered in Tempo2 (TCB 
units). 



MJD" 


P (Aqueye)* 




(s) 


54750.0 


0.033621639166 


54751.0 


0.033621675484 


54752.0 


0.033621711803 



" MJD at the solar system barycenter (Tempo2 mode). 
*" crp= 1. 7 ps (68% statistical error). 



* TCB is a coordinate time referred to the barycenter of the Solar 
System, synchronized with the proper time of a distant observer comov- 
ing with it. The system adopted in Tempol is the barycentric dynamical 
time (TDB), effectively measured in units that differ subtly from the 
conventional SI second (Hobbs et al. 2006). It is as if the time dilation 
effects were not correctly accounted for using TDB units, so that, for ex- 
ample, rotational periods in the TDB system are systematically shorter 
than the TCB ones. 
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Fig. 3. Difference between the optical and radio time of arrival of 
the main peak of the Crab pulsar The optical peak leads the radio 
one by ~ 230 fis (at MJD=54750, epoch of the first observation). 
The (blue) line is the radio-optical drift, which is consistent with 
zero within the errors (see text). 




-100 -50 50 100 

dt (/is) 



Fig. 5. Same as Figure |4] for a simulated signal with superim- 
posed random noise (see text for details). The assumed count rate 
is the average count-rate measured by Aqueye. The bin width is 
3.6;us. The distribution is fit by a gaussian with cr ~ 24 yus. 




Fig. 4. Distribution of the residual of the optical phase of the 
main pulse around the best fit shown in Figure|2] Each point rep- 
resents the number of phase residuals A^, in the /-th bin. The error 
bar on each bin is V^- The three datasets refer to all the obser- 
vations {circles), to observation 4 {squares) and observation 1 
{diamonds), respectively. The bin widths are 1.5/is, 4.1/is and 
7.3fj.s. The (red) solid lines superimposed to each dataset repre- 
sent the best fitting gaussian with cr ~ 32 //s {circles), cr ~ 24- fis 
{squares) and cr ~ 54 fis {diamonds), respectively. 



4.2. Radio-Optical delay 

The time of arrival at the detector of the first pulse of the light 
curve after a certain epoch f is fa„- = ipPmit, where (p is the phase 
defined in Section [3] Figure |3] shows the difference between 
the optical time of arrival of the main peak of the Crab pulsar 
and the radio one determined from the JB radio ephemerides. 
The dispersion measure at around the epoch of the Aqueye ob- 
servations was 56.7842 pc cm"^ (Oct 15, 2008). We find that 



the optical peak leads the radio one. Taking into account the 
uncertainty on the time of arrival quoted in the JB archiv^ 
(~ 60 fis) and the errors from the fit, the time difference is 
~ 230 + 60 fis at MJD=54750, with a drift of ^ 7 fis/day. While 
within the errors the drift is consistent with zero, the difference in 
the arrival times is significant. The optical peak leads the radio 
one, in agreement with what wa s fo und previously b y 'Sanwall 
(ll999h . lOosterbroek etalJ (l2006h and lOosterbroeke t al. (20oC 
the latter obtained using simultaneous optical and radio obser- 
vations. Our value of the radio delay is also consistent with 
the rece nt measurement per formed by our group with Iqueye 
(1 78 us; Colhns et al. 12012b . but larger than the one reported 
in IShearer et al. ( 2003h . The uncertainity on our measurement 
is dominated by the error on the radio ephmerides, and can be 
further reduced in future using simultaneous radio-optical ob- 
servations. 



4.3. Phase noise 

Figure|4]shows the distribution of the phase residuals in Figure|2] 
The distribution is fit with a gaussian and, for all the observa- 
tions, gives a reduced;t'^ ~ 0.9. The gaussian has cr ~ 32 fis, con- 
sistent with the error bar of each measured phase (see Section[3]l. 
This indicates that, with the present accuracy, the phase noise 
of the Crab pulsar observed with Aqueye can be approximated 
with a Gaussian. However, the distribution of phase residuals in 
different observations appears to have different widths. The dis- 
tributions with the smallest (obs. 4) and largest (obs. 1) widths 
are also shown in Figure |4] for comparison. This suggests the 
existence of additional errors in the data chain of some observa- 
tions, possibly induced by signal loss (e. g. clouds or telescope 
tracking errors). 

We compared the observed phase residuals with those ob- 
tained from a synthetic signal generated as a sequence of pulses 
distributed in time with the probability proportional to the pulse 
shape, such as the average count rate and background level 



http ://w w w.j b. man. ac . uk/ pulsar/crab/crab2 . txt 
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are those of the Crab pulsar observed with Aqueye. The sig- 
nal has superimposed random noise and lasts ~ 850 s. Figure |5] 
shows the distribution of the residuals of the simulated signal. 
The distribution is clearly Gaussian but with cr ~ 24fis. This 
is smaller than the total phase noise distribution in Figure H] 
However, as mentioned above, the observations have different 
widths of the distributions. Some of them are comparable (obs. 
4) or marginally larger (obs. 6 and 7) than that inferred from the 
simulation, while some others have more outliers. It is thus pos- 
sible that the outliers are due to errors in the data chain, induced 
by the smaller quality of the dataset. For the best dataset (obs. 4; 
see Figure |4]i, the measured width of the phase residuals distri- 
bution appears to approach the theoretical expectations for phase 
noise induced by pure photon statistics. 

5. Discussion and conclusions 

We observed the Crab pulsar with the photon counting instru- 
ment Aqueye, mounted at the 182cm Copernico telescope in 
Asiago, during the nights of October 10-13, 2008. The counting 
statistics and quality of the data allowed us to monitor the phase 
of the main peak of the Crab pulsar over 2 s long intervals and to 
determine the pulsar rotational period and period derivative with 
great accuracy, using observations covering only a 2 day interval 
in time. The statistical error on the period inferred from a fit of 
the pulsar phase drift is of the order of a few picoseconds. The 
measurements of the period and period derivative agree within 
the statistical error with those inferred from the JB ephemerides. 

We also found that the time of arrival of the optical peak of 
the Crab pulsar leads the radio one (with the latter inferred from 
the JB radio ephemerides) in agre ement with previous findings 
(ISanwal ll999(lShearer et al.l2003h . The actual value of the radio 
delay, ~ 230 yus, is in agreement with the most accurate measure- 
ment p reviously reported in the literature by Oosterbroek et al. 
(l2008f) . Previous measurements performed using fitting func- 
tio ns of the main peak are less pr ecise. In such a case , as pointed 
bv lOosterbroek etal.1 (|2006) and lCollins et all (l2012h . the phase 
of the main peak depends on the chosen fitting function, intro- 
ducing a systematic diff'erence caused by the intrinsic asymmet- 
ric shape of the peak. Assuming a Gaussian or Lorentzian fit, 
different values of the phase of the peak are obtained for differ- 
ent fitting ranges around it. Using a 0. 1, 0.05, 0.02 phase interval 
to the right and to the left of the peak, the difference in position 
is 247 yus, 230 /is, 134 fj.s for the Gaussian fit and 197 yus, 197 //s, 
1 14 yUS for the Lorentzian fit, respectively. Furthermore, it is im- 
portant to note that the correct estimate of the radio-optical time 
delay can be obtained only from simultaneous optical and radio 
observations using si milar procedures fo r the barycenterization 
and the analysis (e.g.lCoUins et al.|| 2012l) . 

As pointed out by Oosterbro ek et al.l (|2006), a time delay of 
the radio peak of ~ 230 fis could have two different interpreta- 
tions. The emission region of the optical radiation is: (a) higher 
in the magnetosphere (~ 70 km) than the radio emission, (b) lo- 
cated at a different angle (~ 2.5") with respect to the radio one. 

We also studied the phase noise distribution of the Crab pul- 
sar observed with Aqueye and found that, with the present ac- 
curacy, it can be modelled as a Gaussian. The width of the dis- 
tribution of the entire dataset is slightly larger than that induced 
by pure photon random noise from a synthetic signal having the 
same pulse shape, average count rate and background level of 
the Crab pulsar observed with Aqueye. However, residual sys- 
tematic errors in the data chain may be present in a subset of 
observations that tend to broaden the distribution. In fact, the 
distribution of the best batch of data, i.e. that with the smaller 



width, is consistent with that induced by photon statistics. While 
the observed broadening seems to be caused mostly by the de- 
terioration in the quality of some observations, the existence of 
a smaller source of phase noise, possibly related to the intrinsic 
pulsar mechanism, cannot be ruled out at present and needs to 
be carefully investigated with future observations. 
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Appendix A: Calculating the phase of the main peak 

We start folding the light curve over a fixed reference period 
Pi„i,, corresponding to the period of the pulsar some time during 
the complete observing run. 

The phase drift of the main peak i// (see Section[3]l is obtained 
by finding the maximum of the coiTelation function K between 
the actual pulsar signal P{(p(t)) and a template S{t/Pi„it - x) with 
a time dependence closely resembling the actual pulsar signal, 
but with a varying phase x. The function S (x) is a periodic func- 
tion with period 1 , obtained as a smoothed version of the whole 
dataset (using 3504 data points in phase): 



where At is the averaging time (usually 2 seconds) and t\ + At/2 
the mid point of the averaging interval. After taking the deriva- 
tive of K with respect to x, we obtain: 



The value of is obtained taking K'(ti,x) - and considering 
only the zero with the steepest crossing (see Fig. lA.ll ). 

We select the steepest zero crossing because we are look- 
ing for the maximum which corresponds to the sharpest peak of 
K. Let f(x) be a differentiable function with many extrema xi, 

X2, Then, in the vicinity of each peak it can be expanded as 

fix) - f{xk) + {\ I2)f"{xk){x - Xk)^ + The function has maxi- 
mum at Xk if f"(xk) < and a minimum if f"{xk) > 0. f'{x) has 
a zero at Xk and the slope of f'(xk) is f"(xk). The absolute value 
of f"{xk) - f{xit)/w^, where w is the width of the peak of the 
function at half maximum. Since it is unlikely that a correlation 
function would have a small very narrow peak that would make 
I/" I larger than the main peak, one concludes that the highest 
peak also has the largest value of |/"| and hence the steepest |/'|. 
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Fig. A.l. Example of the behaviour of K (f i , x). 
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